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It is argued that in relativistic heavy ion colhsions, due to limited size of the formed matter, 
the rehable criterion of critical point is finite-size scaling, rather than non-monotonous behavior of 
observable. Ifow to locate critical point by finite-size scaling is proposed. The data of pt correlation 
from RHIC/STAR are analyzed. Critical points are likely observed around y/s = 62 and 200 
GeV. They could be, respectively, the transition of deconfinement and chiral symmetry restoration 
predicted by lattice-QCD. Further confirmation with other observable and energies is suggested. 
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Lattice-QCD simulations have shown that the tran- 
sition of deconfinement in quantum chromodynamics 
(QCD) at vanishing baryon chemical potential /^b is 
crossover There has been much speculation that the 
crossover becomes a true first-order phase transition for 
larger values of ^b- This suggests that the QCD phase 
diagram can exhibit a critical endpoint where the line 
of first order transition matches that of second order or 
analytical crossover 0. 

Chiral symmetry restoration is another QCD origi- 
nated phase transition. It has been shown that the tran- 
sition for yUB = is crossover So there could also 
be a chiral critical endpoint in phase diagram. But it 
is unclear if the critical temperature of chiral symmetry 
restoration is above [4], or equal to ^5], or below [gj] that 
of the deconfinement. 

Locating the critical endpoints of QCD phase transi- 
tions by lattice calculation is still a formidable challenge. 
But if the critical endpoint is in the region accessible to 
current relativistic heavy ion collisions, it should be dis- 
covered experimentally. 

Most of the current signatures for finding the critical 
point are focused on the anomalous, or non-monotonous, 
behavior of the observable at various incident energies 7] . 
The argument is that in infinite system, the correlation 
length ^ diverges when approaching the critical point. 
The contribution of this singularity to the observable is 
supposed to be proportional to However, the data 
from RHIC and SPS in more than a decade accumula- 
tion show no sign of anomalous behavior as a function of 

75 i. 

In relativistic heavy ion collisions, two nuclei move 
with relativistic velocity and collide as two contracted 
pancakes. More central collision makes overlapped area 
larger. It is just because the large number of strongly in- 
teracting nucleons in more central nuclear collisions make 
the transition between hadron and quark-gluon plasma 
possible. The centrality (or the system size) dependence 
of the observable is noticeable 0]. 

Due to the finite size of system, no divergence can be 
practically observed at critical point. The physical quan- 



tities, which are divergent in infinite system, become fi- 
nite and have a maximum, i.e., so called non- monotonous 
behavior. However, the position of the maximum changes 
with system size and deviates from the true critical point. 

The appearance of non-monotonous behavior is not 
always associated with critical point. Taking one- 
dimensional Ising model as an example, there is no crit- 
ical point in this model, but its specific heat in a finite 
system has non-monotonous behavior. 

Moreover, the absence of non-monotonous behavior 
does not mean no critical point. The physical quantities 
like order parameter, which are finite in infinite system, 
have a monotonous behavior near critical point in a finite 
system (loj . Therefore, non- monotonous behavior is not 
a reliable criterion for the critical point of finite system. 

An effective identification of critical point of finite sys- 
tem is the finite-size scaling, which was proposed from 
phenomenological [ll| and renormalization-group [l^ 
theories, and was approved by the Monte Carlo results 
of finite systems in different universal classes [l^ . 

In this letter, we first propose how to locate critical 
point by finite-size scaling. Then the data of pt cor- 
relation at 6 centralities and 4 incident energies from 
RHIC/STAR are analyzed. The behavior of critical point 
is likely observed around -/s — 62 and 200 GeV. Finally, 
we suggest how to confirm the findings and precisely lo- 
cate the critical point in coming experimental study. 

The main points of finite-size scaling can be described 
as the following. An observable Q of finite system is a 
function of temperature T and system size L. When L is 
much larger than the microscopic length scale and T is 
in the vicinity of critical point Tc, the observable Q(T, L) 
can be written in a finite-size scaling form [ill ll2L \13^ , 

QiT^^^L^'^FQitL^'"). (1) 

t ~ (T — Tc)/Tc is the reduced temperature and A is 
the critical exponent of the observable, v is the critical 
exponent of the correlation length ^ = ^o^"''- 

Finite-size scaling not only characterizes the scaling 
behavior of thermodynamic quantities of finite system 
near critical point, but also provides criterion for locating 



2 



the critical point. At critical point T = Tc, the finite-size 
scaling function Fq in Eq. ([T]) becomes 

FQ(0) = Q(^e,L)L-^/^ (2) 

which is constant and independent of system size L. In 
the plot of Q{T,L)L-^/'' vs T, the critical point [Tc, 
Fq{0)] is a fixed point, where all curves of different sys- 
tem sizes converges to. Reversely, the appearance of fixed 
point indicates the existence of a critical point. 

If the critical exponent A = 0, like Binder cumulant 
ratio the fixed point can be obtained directly from 
the temperature dependence of this observable at differ- 
ent system sizes. This is why Binder cumulant ratio has 
been used very widely in determining critical point of 
finite-size system. 

If the critical exponent A 7^ and is unknown, the 
fixed point can be found by investigating the tempera- 
ture dependence of Q{T, L)L^"^ at different system sizes. 
When a fixed point is observed at a certain parameter 
ao, it indicates the existence of a critical point and the 
parameter ao is related to the ratio of critical exponents, 
i.e., X/u = flQ. 

The critical point can also be found directly from the 
system size dependence of the observable. Taking loga- 
rithm in the both sides of Eq. ([T]), it becomes 

\nQ{T,L) = X/vlnL + lnFQitL^/"). (3) 

At critical point t = 0, the second term of Eq. ^ be- 
comes a constant and lnQ{Tc, L) becomes a straight line 
with respect to InL. If system is away from the critical 
point, the second term of Eq. ([3]) is no longer a con- 
stant. It gives an additional size dependent contribution 
to the observable and makes \nQ{T, L) deviate from the 
straight line with respect to In L. 

It is found recently that the finite-size scaling holds not 
only for thermodynamic quantities like order-parameter, 
susceptibility, and so on, but also for various cluster 
sizes [1^ and their fluctuations [l^. Therefore, the finite- 
size scaling of various critical related observable could be 
used to identify critical point and its critical exponents. 

In relativistic heavy ion collision, correlation and fiuc- 
tuation of final state particles is regarded as critical re- 
lated observable Although much attention have 
been drawn in measuring them, but influence of system 
size has been neglected. The available data for system 
size study is very few. The pt correlation at Au-I- Au 
collisions from RHIC/STAR Q is the only data which 
can be used for the analysis, where the centrality de- 
pendence of Pt correlation at 4 incident energies are well 
presented 0] . But the errors of the data at y/s = 20 GeV 
are much larger than that at other collision energies. The 
Pt correlation is defined as 

^ 7v„ E E {pt,t - {pt)){pt.j - (pt)) 



Nc is the number of event, pt,i is the transverse- 
momentum of the ith particle in each event, and Nk is 
the number of particles in the kth event. (...) is the 
average over event sample. 

Collision energy is the controllable condition. Here we 
let it play the role of temperature in the analysis of finite- 
size scaling. The size of the formed matter is mainly lim- 
ited by the size of overlapping transverse region, which is 
proportional to the number of participant nucleons and 
is quantified as centrality. So the initial mean size of the 
formed matter can be approximately estimated by the 
square root of participants, -^/iVpart • We choose dimen- 
sionless (or relative) size, 

L = V^V;^/v/27Va, (5) 

as scaled mean size of initial system, where Na is the 
number of nucleons of incident nucleus. The system size 
at transition should be a monotonically increasing func- 
tion of L. The position of critical point is insensitive to 
the concrete form of this function, but only the critical 
exponents changes with it. As the first step, the system 
size at transition is assumed to be proportional to L. 

In the case of a few critical points, the finite-size scaling 
of pt correlation in the vicinity of each critical collision 
energy y/sZi{i — 1, 2, ...) can be written as 

P(yi,L)=L^'/'''Fp,,[e,Li/''']. (6) 

= (a/s— y/Sc.i)/ \/Sc,i is the reduced collision energy at 
ith critical point, which is unknown in priori. Ai is the 
ith critical exponent of pt correlation. In the following, 
we demonstrate how to locate the critical point by the 
data of Pt correlation from RHIC/STAR. 

Firstly, we change the centrality dependence of pt cor- 
relation at different collision energies in Ref. [9| to the 
collision energy dependence at different sizes (or central- 
ities). The results are shown in Fig. 1(a). Since in the 
most peripheral collisions, the size of the formed matter is 
too small to be inside the asymptotic region of finite-size 
scaling, we choose six centralities at mid-central and cen- 
tral collisions to do the analysis. The sizes correspond- 
ing to the 6 centralities are indicated in the legend of 
Fig. 1(a). It is clear that at a given collision energy, the 
correlation strength increases with the decrease of system 
size. The influence of finite size is obvious. 

If critical collision energy of QCD phase transition is 
in the range of incident energy at RHIC, the behav- 
ior of fixed point should be observable. So we multiply 
P(-y/s, L) by a size factor with different a to see how 
it changes with the system size L. Varying —a from small 
to large, it is interesting to see that at collision energy 
-y/s = 62 GeV, all points of different sizes move firstly 
toward each other, then well converge at — ao.i — 2.09, 
and finally move again apart from each other. The cor- 
responding steps and typical a values are presented in 
Fig. 1(b), (c), and (d) respectively, where the errors in 
each sub-figures come from the measure of P{\/^, L) only, 
and the errors of A^part are not included. 
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FIG. 1: (a) The energy dependence of pt correlations at different sizes L (or centralities) . Data come from RHIC/STAR [1). 
(b), (c) and (d) are pt correlation multiplied by the factor, I/"", with —a = 1.0, 2.09 and 4, respectively 



TABLE I: Parameters of parabola fits. 



Vs(GeV) 


20 


62 


130 


200 


|C2| 


1.86± 0.93 


0.6 ± 0.09 


1.56± 0.41 


0.77± 0.1 




3.9±0.89 


2.59± 0.09 


3.43± 0.41 


2.74±0.1 



At y/s = 200 GeV, the points of different sizes show 
the same behavior and best converge at — ao,2 = 2.08. 
While in the whole process, the points of different sizes 
at energies -^/s = 20 (or 130) GeV never move close to 
each other as those at y/s — 62 (or 200) GeV do. So 
there are likely two fixed points around y/s = 62 and 200 
GeV. 

In order to confirm the position of fixed points, we 
study the In L dependence of In P{\/s, L) for four incident 
energies, respectively. A parabola fit, C2(lnL)^+ci lnL + 
Co, is used at each collision energy. The better straight- 
line behavior results in smaller \c2\ and larger ratio of 
|ci/c2|. The fit parameters, C2 and ci, for 4 collision 
energies are listed in Tab. 1. It shows that the better 
straight-line behavior happen to be at y/s = 62 and 200 
GeV, which arc the same collision energies of fixed points 
found above. The data at these two energies can be well 
fitted, respectively, by the straight lines with slopes ao,i 
and ao,2 obtained above by the fixed points. The results 
are shown in Fig. 2(a). While, the data at y/s — 20 
and 130 GeV are better fitted by parabola as shown in 
Fig. 2(b). 

The same analysis has also been applied to the pt corre- 
lation normalized by the average pt over the whole sam- 
ple i. The analysis for normalized pt correlation at 
-y/s = 62 and 200 GeV show exactly the same behavior 
of fixed points and straight lines as what pt correlation 
demonstrates above. The critical exponents of normal- 
ized pt correlation (about 1.1) are smaller than that of 
Pt correlation. 

So the critical collision energies are most probably 
around ^/s = 62 and 200 GeV, rather than near y/s — 20 
and 130 GeV. The same analysis for other critical related 
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FIG. 2: Double-log plots of pt correlation with respect to 
size, (a): straight-line fits with slopes ao,i and ao,2 obtained 
by fixed points, and (b); parabola fits. 



observable, such as the fluctuations of mean pt per event, 
the moments of multiplicity, the ratio of K to tt, and so 
on, will be greatly helpful in confirming the observed re- 
sults. Therefore, the incident energy and centrality de- 
pendence of those observable are called for. 

If there were additional collisions around ^/s — 62 and 
200 GeV, we could determine the finite-size scaling func- 
tion defined in Eq. ([6]). This is impossible at present 
since there are only two collision energies in addition to 
the critical ones, and they could be outside of the asymp- 
totic region where finite-size scaling holds. 

The findings of the two critical points may imply that 
deconfinement and chiral symmetry restoration occur at 
different temperatures. Which one is at the lower or 
higher temperature (energy) has to be confirmed finally 
from theoretical calculation. Two critical collision ener- 
gies, y/s — 62 and 200 GeV, are both within the range 
estimated by lattice calculation (isj . 

The similar ratios of critical exponents at two criti- 
cal points is consistent with current theoretical estima- 
tion, which shows that all critical exponents of the de- 
confinement transition, in the same university as the 3- 
dimensional Ising model [3, are very close to that of 
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chiral symmetry restoration, in the same university as 
the 3-dimensional 0(4) model with spin symmetry [20|. 

To the summary, we argue in this letter that finite- 
size effects of the formed matter in relativistic heavy ion 
collisions is not negligible. The finite-size scaling, rather 
than non-monotonous behavior of observable is a reliable 
criterion of the existence of critical point. Then we pro- 
pose how to locate critical point by finite-size scaling. As 
an application, we analyze the data of pt correlation and 
its normalized one at 6 centralities and 4 incident ener- 
gies from RHIC/STAR. Two fixed points, and therefore 
two critical points, are likely observed around ^/s = 62 
and 200GeV. They could be, respectively, related to the 
transition of deconfinement and chiral symmetry restora- 
tion predicted by latticc-QCD. The ratios of critical ex- 



ponents at these two critical points are similar, in con- 
sistence with current theoretical estimation. 

The confirmation of this observation requires the ef- 
forts from both theoretical and experimental sides. From 
experimental side, it is proposed to get more and bet- 
ter data on other critical related observable at current 
collision energies, and a few additional collisions around 
y/s — 62 and 200 GeV. Then we can more precisely de- 
termine the critical cndpoints and critical exponents. 
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